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Abstract 

In the framework of the perturbative QCD approach, we study the charmless pure weak anni- 
hilation B- -> K~K° decay and find that the branching ratio BK(B C ^KK) ~ 0(1(T 7 ). This 
prediction is so tiny that the B c — > KK decay might be unmeasurable at the Large Hadron Collider. 
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I. INTRODUCTION 



The study of the decays of B mesons is important and interesting for the determination 

of the flavor parameters of the Standard Model (SM), the exploration of CP violation, 

the search of new physics beyond SM, etc. In recent years, theoretical studies of B u d 

mesons have been investigated widely in the literatures. They are tested and supported by 

the experimental data collected by the detectors at the e + e~ colliders, such as the CLEO, 

Babar, and Belle. With the bright hope arising from the startup of the CERN Large Hadron 

Collider (LHC) [l| , the heavier B s and B c mesons could be produced abundantly and studied 

in detail at the hadron colliders. It is estimated that one could expect around 5 x 10 10 "self- 

f] 

■lagging" B c events per year at the LHC [2| due to a relatively large production cross section 
3] plus the huge luminosity C = 10 34 cm~ 2 s _1 [4]. There seems to exist a real possibility to 
study not only some B c rare decays, but also CP violation and polarization asymmetries. 
The study of the B c mesons will highlight the advantages of B physics. 

The B c mesons are the "double heavy-flavored" binding systems and share many features 
with the heavy quarkonia. The first observation of the B c mesons at the Tevatron jlj] 
provokes the physicist's particular interest in them. Many studies and investigation of the 
properties of the B c mesons have been made, and will be further scrutinized by the LHC 
experiments. Because the B c mesons lie below the BD threshold (here we only discuss the 
lightest I^Sq ground state pseudoscalar B c mesons, excluding their excited states) and carry 
flavors, they cannot annihilate into gluon and/or photon so are stable for the strong and/or 
electromagnetic interaction. Because of the flavor quantum numbers B = — C = ±1, the 
B c mesons can decay through the weak interaction only. The B c mesons have more decay 
modes than the B u ^ s mesons due to several reasons. One is that many decay modes, such 
as B c — > B u ^ s + X, are only accessible by B c mesons because of their sufficiently large 
masses. Another is that the B c mesons carry open flavors, so either b or c quarks can decay 
individually. The potential decays of the B c mesons permit us to over-constrain quantities 
determined by the B u ^ s meson decays. 

The decays of B c mesons can be divided into three classes: (1) the 6-quark decay (i.e. 
b — >■ q u , where the up-type quark q u = u, c) accompanied with the spectator c-quark, 
(2) the c-quark decay (i.e. c — >• q D , where the down-type quark q D = d, s) accompanied 
with the spectator 6-quark, and (3) the annihilation channel (i.e. B~ — > i^ue, Q D %i where 



2 



the lepton i~ = e~ , pT , t~). Among the multitudinous B c decays, the weak annihilation 
channels are expected to take ~ 10% shares according to the estimates in 2|, |6( for which a 
major part comes from the tree weak annihilation process B~ — > sc which is not helicity- 
suppressed because of the large charm quark mass and produces a large weak annihilation 
branching ratio with charm in the final state, while the charmless pure weak annihilation 
decay B c — > KK is helicity-suppressed like the Bd — > KK decay and would have a very 
small branching ratio. It is highly expected that the LHC experiments might shed light on 
a better understanding of weak annihilation processes for B c mesons. 

In recent years, several attractive methods have been proposed to study the nonleptonic 
B decays, such as the QCD factorization (QCDF) 7|, perturbative QCD method (pQCD) 
js-lOj, soft and collinear effective theory 11, 12|, etc. Here, we would like to investigate the 
charmless pure weak annihilation B c — > KK decay with the pQCD approach due to several 
reasons. (1) One reason is that fits of nonleptonic charmless decays B u d — > PP, PV without 
taking into account weak annihilation contributions are generally of poor quality [l^l (here 
P and V denote the lightest ground pseudoscalar and vector mesons, respectively). Our 
present understanding of the weak annihilation contributions remains limited and unclear. 
So the pure weak annihilation processes, such as B c — > KK decays, are interesting and 
worthy of study, which will certainly help us to improve our understanding of the weak 
annihilation contributions. (2) Another is that due to both kinematic improvement from 
the large phase spaces and dynamic enhancement of the CKM factor |Kf,V^* d |, the B c — > KK 
decay is expected to have a large branching ratio among two-body nonleptonic charmless 
W- annihilation B c — > PP processes. In addition to the absence of penguin operators for 
the tree annihilation process B c — > KK, final state interactions arising from soft gluon 
exchanges are expected to be extremely small because of the large momenta of the final K 
mesons. Therefore a relatively accurate estimation of annihilation contributions could be 
obtained effectively from the charmless B c — > KK decay. (3) Still another is that Ref.fl^j 
obtains a very large B c — > KK branching ratio, about 1.6%, at 4 orders of magnitude 
bigger than the estimate O(10~ 6 ) of Ref. 15], but this estimate is not valid because Ref. [jjj] 
in their calculation incorrectly uses the measured penguin-dominated B^^tn^K branching 
ratio while the decay B c — > KK is a pure tree weak annihilation and should be related to 
Bd — > KK. In addition, the branching ratio of the charmless decay B c — > KK is estimated 
to be O(10~ 8 ) with the QCDF approach 15j. Recently, this charmless decay is also studied 
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with the pQCD approach and its branching ratio is O(10 7 ) with the off-mass-shell final 



states 



161 ]. which is the same order of magnitude as ours obtained in this paper with the 



on-shell final states. 

This paper is organized as follows : In Section HI1 we will discuss the theoretical framework 
and give the decay amplitudes for B c — > KK with the perturbative QCD approach. In our 
calculation, we shall ignore the final state interactions because the final states have very 
large momenta and move far away before soft gluon exchange. Section HTT1 is devoted to the 
numerical result of the branching ratio. Finally, we summarize in Section IIVI 



II. THEORETICAL FRAMEWORK AND THE DECAY AMPLITUDES 



A. The effective Hamiltonian 



Using the Operator Product Expansion approach and renormalization group (RG) equa- 
tion, the low-energy effective Hamiltonian for B c — > KK decay can be written as 

«.// = |w{CiMQi + CMQ 2 ) + H.c, (1) 

where Gf is the Fermi coupling constant for electroweak interactions. V c bV* d is the CKM fac- 
tor accounting for the strengths of the nonleptonic B c decays. Cj(/i) are Wilson coefficients 
at the renormalization scale /i which have been evaluated to the next-to-leading order with 
the perturbation theory. The local tree operators are process dependent. Their expressions 
are defined as 

where a, ft are SU(3) color indices. The most difficult problem in theoretical calculation 
of nonleptonic charmless decay B c — > KK is how to evaluate the hadronic matrix elements 
(KK\Qi t2 \B c ) properly and accurately. 



B. Hadronic matrix elements 



For convenience, the kinematics variables are described in the terms of the light cone 
coordinate. The momenta of the valence quarks and hadrons in the rest frame of the B c 
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meson are defined by 

Pb~ = Pi = ^(!> !> 0j-)> fc i = x iPi + (°> °> 

Pk- = P2 = ^(1, 0, 0i), fc 2 = x 2 p 2 + (0, 0, A? 2 ±), n 2 = (1, 0, 0), 

Pk<> = P3 = ^(0, 1, Ox), h = x 3 p 3 + (0, 0, k 3± ), n 3 = (0, 1, 0), 

where n2-n 3 = 1. The null vectors n 2 and n 3 are the plus and minus directions, respectively. 
k\ is the momentum of c quark in the B c meson. k 2 and k 3 are the momenta of the light 
non-strange quark in the K~ and K° mesons, respectively. kn_ denotes the transverse 
momentum. X{ denotes the longitudinal momentum fraction of the valence quark. 

The calculation of the hadronic matrix elements is difficult due to the nonperturbative 
effects arising from the strong interactions. Phenomenologically, using the Brodsky-Lepage 



approach 



17J, a modified perturbative QCD formalism has been proposed recently under 



the kr factorization framework [8l4l0j|. Taking into account the transverse momentum of 
the valence quarks in the hadrons, the Sudakov factors are introduced to modify the end- 
point behavior of the hadronic matrix elements. The amplitudes are factorized into three 
convolution parts : the "harder" functions, the heavy quark decay subamplitudes, and the 
nonperturbative meson wave functions, which are characterized by the W ± boson mass m\y, 
the typical scale t of the decay processes, and the hadronic scale Aqqd, respectively. The 
pQCD approach has been extensively applied to study semileptonic and nonleptonic B de- 
cays with phenomenological results. More information about pQCD approach can be found 
in [8l4l0l|. The final decay amplitudes can be expressed as 

A(B-^K~K Q ) oc C{t)®H{t)®$ B - {xi, h)^ K - (x 2 , h)®$ K o(x 3 , b 3 ), (3) 

where the Wilson coefficient C(t) is calculated in perturbative theory at the scale of mw 
and evolved down to the typical scale t using the RG equations. <g> denotes the convolution 
over parton kinematic variables. H(t) is the hard-scattering subamplitude which is domi- 
nated by hard gluon exchange and can be factorized. The universal wave functions $(2,6) 
absorb nonperturbative long-distance dynamics, which can be extracted from experiments 
or constrained by lattice calculation and QCD sum rules, b is the conjugate variable of the 
transverse momentum of the valence quark of the meson. According to the arguments in 



lOj , the amplitude of Eq. ([3]) is free from the renormalization scale dependence. 
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C. Bilinear operator matrix elements 



Within the pQCD framework, the long-distance hadronic information is contained by the 
the so-called light-cone distribution amplitudes (LCDAs) which are defined from hadron-to- 
vacuum matrix elements of nonlocal bilinear operators. Although LCDAs are not calculable 
in QCD perturbation theory, some of their properties are well understood for both light and 
heavy mesons. For example, the LCDAs for the K meson including higher-twist contribu- 



tions are systematically presented in (18J. In our calculation, we only consider two-particle 
(valence quarks) twist-2 and twist-3 LCDAs for K mesons, and neglect contributions from 
higher Fock states. The LCDAs for K mesons are written as 

(if(p)|s Q (0)g^)|0}=-^ (4) 

where N c is the color number. The parameter fix is the chiral factor = m 2 K /(m s + m q ). 
The null vector n + and n_ are parallel to p and z, respectively. The expressions of the 
twist-2 LCDAs (p a K and the twist-3 LCDAs (p p K , <$ K are collected in Appendix IA1 

Unlike the ir and K mesons, our knowledge of the LCDAs for B c mesons has been 



relatively poor until recently (for a recent view, see 19j). but we know that the B c mesons 
are composed of heavy valence quark both b and c. Given m# c w + m c , the B c mesons 
can be described approximately by nonrelativistic dynamics. In this paper, we will take 

(0\c a (z)b^0)\B:(p 1 )) = ^Id^e-^-f^+m^Ts^^)]^, (5) 

where Jb c is the decay constant of the B c meson. As the arguments in Jj|, this simplest 
form, <Pb c { x ) = ~ m c/^B c )i is the two-particle nonrelativistic LCDAs at the tree level 
where both heavy valence quarks just share the total momentum of the B c mesons according 
to their masses. For a rough estimation of the branching ratio for B c — y KK decay, we will 
take the simplest form as an approximation, and neglect the relativistic corrections and 
contributions from higher Fock states. 

D. The decay amplitudes 

The B c — y KK decay is the pure annihilation process. According to the effective Hamil- 
tonian Eq.Q, the lowest order Feynman diagrams are shown in FIGJTJ where (a) and (b) 
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are nonfactorizable topologies, (c) and (d) are factorizable topologies. After a straightfor- 
ward calculation using the modified perturbative QCD formalism Eq.(j3]), we find that the 
contributions of factorizable topologies are zero, which is a result of exact isospin symmetry. 
The decay amplitude comes only from the nonfactorizable topologies, and can be written as 

= ~ i TkIT - V cbVud dxidx 2 dx 3 / &id&i / b 2 db 2 (f) Bc (x 1 ) 

y2N c Jo Jo Jo 

x {a s (t a )C 2 (t a )E(t a )H(A, OC, h, b 2 ) \(f>K- ( l ) K ( r b + x l — x 3) 
+ r K -r K o(jf K -4F KO (Ar b + 2x\ - x 2 - x 3 ) 
+ r K -r K o (4>k-4>k° + ^k-^ko) ( x 3 - x 2 ) 
+ r K -r K °<t>K-(f>Ka {2xt -x 2 - x 3 ) 

+ a s (t b )C 2 (t b )E(t b )H(A, 0, b\, b 2 ) 



F K -(j) a K ° O2 -xi- r c ) 



+ r K -r K o(f) p K _(j) p KO (x 2 + x 3 - 2xi - 4r c ) 
+ r K -r K o (f K -()f K0 + ( f K -4 >t Ko) (^3 - x 2 ) 

+ r K -r K o(j) t K -(f) t KO (x 2 + x 3 -2x 1 )}} (6) 

j j 62=03 

where the CKM matrix elements V cb V* d = AX 2 (1 - \ 2 /2 - A 4 /8) + C(A 8 ) with the phe- 
nomeno logical Wolfenstein parameterization. r b = m b /mB c and r c = m c /mB c are the ratios 
of the mass of b and c quark to the mass of B c mesons, respectively, rx = ^k/^b c = 
m 2 K /[mBX m s Cf = 4/3 is the SU(3) color factor. t a ( b ) is the characteristic scale. 

A is the virtualities of internal gluons, which is a timelike variable for the pure annihilation 
B c — > KK decay concerned, a and /3 are the virtualities of internal quarks. E and H are 
the Sudakov factor and the hard kernel functions, respectively. Their expressions are listed 
in Appendix iBl 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The branching ratio in the B c meson rest frame can be written as: 



BK(B C ^KK) 



tb c V 
8tt ra 2 B 



■\A{B C -*KK)\' 



(7) 



where p is the center-of-mass momentum of K mesons. The lifetime and mass of the B, 



meson are m Bc = 6.276 ± 0.004 GeV and r Bc = 0.46±0.07 ps 



22j, respectively. Other input 



parameters are 

m c = 1.2718;2I GeV [22], A = 0.2257±g;8g?g [22], / Bc = 489±4 MeV [23], 

m fc = 4.201&07 GeV [22], A = 0.814lg;8^ [22], /«■ = 159.8±1.4±0.44 MeV [24]. 

If not specified explicitly, we shall take their central values as the default input. The nu- 
merical result of the branching ratio is 

Bn(B;^K-K Mimt f 7 (m b )t° ^ 

where the errors come from the uncertainties of quark masses m& and m c , the CKM factor 
V c bV* d , and the decay constants fs c and fx- The largest error arises from the parameter of 
m&, which can reach 40%. The errors arising from both the CKM factor and the decay con- 
stants are relatively small. Of course, there are some other uncertainties not considered here, 
such as the radiative corrections to the LCDAs of B c mesons, the final states interactions, 
etc. So the results might just be an estimation of the pQCD approach. 

Our estimation of the branching ratio B1Z{B C -^KK) is slightly different with the result 
in [if]], although they are calculated with the same pQCD approach resulting in the same 
order of magnitude O(10~ 7 ). Besides the input parameters, the reasons may be (1) whether 
the final states are on-mass-shell or not, and (2) whether the contributions of factorizable 
topologies are zero or not c . With appropriate input parameters, the results in [3] and ours 
are in agreement with each other within an error range. 

As the arguments in [ijj], the inconsistencies among various estimations of the branching 
ratio BTZ(B C -^KK), such as O(10~ 6 ) based on annihilation by using the relations among 
;he charmless weak annihilation B c decay channels relying on the SU(3) flavor symmetry 
lf| , O(l0~ 7 ) (or O(10~ 8 ) [ljj]) based on perturbative one-gluon exchange with the pQCD 
(or QCDF) approach, arise from conceptually different methods. Anyway, for weak annihi- 
lation to light quarks in the final state, the tree annihilation B~ — > du process is helicity 
suppressed because of small light quark masses, so that gluon emission either from the ini- 
tial or final state must occur in this annihilation and the decay amplitude is then 0(a s ) 



as given in pQCD. Both the estimations in 



15 



161 ] and our result are in accordance with 



Because of almost equal masses of the final states, the same two-particle LCDAs for the charged and 
neutral K mesons are taken in our calculation. With this approximation, a similar conclusion, that the 
contributions of factorizable topologies cancel each other because of the isospin symmetry, can be found 
in B s — > 7T7T decays with the pQCD approach jgsj. 
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an intuitive expectation for nonleptonic charmless ^-annihilation of heavy meson decays 
which are usually suppressed. There are some additional factors for the tiny estimation of 
B1Z{B C -^KK). One the is that although the B C -^KK decay is a tree weak annihilation 
process, its amplitude is color suppressed and associated with C<ijN c . Another is that there 
is a large destructive interference between the nonfactorizable topologies due to the near 
equal final state particle masses. This can be clearly found in Eq.([6]). The numerical results 
also confirm the cancellation between the nonfactorizable topologies, and give the strong 
phases ~ —31° and ~ +127° for FIG{1] (a) and (b), respectively. 

If the pQCD prediction is right, then there should be some 10 3 events for B c -^rKK decay 
per year at the LHC Considering the detection efficiency and selection efficiency, there 
would be just a few events per year. The signal of the pure weak annihilation B c -^rKK decay 
would be very tiny at the LHC. As B nonleptonic charmless decays, the charmless pure weak 
annihilation is expected to be small in B c nonleptonic decays, so the LHC measurement could 
confirm our understanding of the annihilation terms in weak decays based on perturbative 
QCD. 

IV. SUMMARY 

In this paper, we study the B c — > KK decay with the pQCD approach, which would 
call for another reassessment of the weak annihilation processes and might provide some 
valuable hints of our understanding on perturbative QCD and long-distance contributions. 
It is found that the contributions of factorizable annihilation topologies are zero, and that 
there is a large cancellation between the nonfactorizable topologies, which result in the 
branching ratio B1Z(B C -+KK) ~ O(10~ 7 ). The branching ratio with the pQCD approach 
is so tiny that the B c — > KK decay might not be measured at the LHC experiments. 
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Appendix A: Distribution amplitude of the K meson 

The expression of the LCDAs of the K meson induing higher-twist contributions can 
be found in 



181 . In our calculation, the twist- 2 distribution amplitude S a K and the twist-3 

n 

distribution amplitude 4> P K and (\> l K are [2CJ 

twist-2 <j> a K {x) = — ^=6xx{l + 0.17C^ /2 (t) + 0.11BC 2 /2 (t)} (Al) 

twist-3 <jf K {x) = ^={l + 0.24C 2 1/2 (t) - 0.12Cl /2 (t)}, (A2) 

twist-3 0Sc(x) = ^^{c} /2 (t) + 0.35Cj /2 (t)} (A3) 

where the decay constant fx = 160 MeV. t = x — x — 2x — 1. The Gegenbauer polynomials 
are 

C't /2 (z) = 3z, C 3 2 /2 (z) = ^(5z 2 -l), 

Cl /2 (z) = z, C 1 2 /2 (z) = -(3z 2 -l), 



Cl ,2 {z) = -(5z 3 - 3z), C\' 2 {z) = -(35z 4 - 30z 2 + 3) 
2 8 

Appendix B: Some parameters and formulas 

The expression of the Sudakov factors E is 

E(t) = exp (S Bc (t) - S K -{t) - S K o(t)) (Bl) 



where 



S Bc (t) = s(x 1 p+,b l ) + 2[ t % q (B2) 

J l/bi /I 

S K -(t) = s(x 2 pt, b 2 ) + s(x 2 p+, b 2 ) + 2[ t (B3) 

J 1/62 /i 

/t 
lq (B4) 
l/£>3 M 



1/63 

The anomalous dimension of the quark is 7 g = —a s /ir. The explicit expression of s(Q,b) 



can be found in 



2l|. 



The hard kernel function H is defined as follows 



H(A,Z,b 1 ,b 2 ) = i [ e(b 1 -b 2 ) l -^Hj ) 1) (VAb 1 )J (VAb 2 ) + (b 1 ^b 2 ) 

x {e(Z)K (VZh) + ei-^jH^iy/lZlbx)} (B5) 
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where the hard scales are 

A = m 2 (l -x 2 )(l -x 3 ) (B6) 

—a = m 2 (xi — X2){x\ — £3) — ml (B7) 

-P = m 2 (x l + x 2 - l)(xt + x 3 - 1) - m 2 c (B8) 

t a = max (VA, l/&i, l/fc) (B9) 

t 6 = max (v'A, ^0, l/&i, l/b 2 ) (BIO) 
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FIG. 1: The lower order feynman diagrams contributing to the B c — > KK decay, with (a) and (b) 
for nonfactorizable annihilation, (c) and (d) for factorizable annihilation. 
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